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Kinetic Vlasov Simulations of collisionless magnetic Reconnection 

H. Schmitz and R. Grauer 
Theoretische Physik I, Ruhr-Universitat, 44780 Bochum, Germany 

A fully kinetic Vlasov simulation of the Geospace Environment Modeling (GEM) Magnetic Re- 
connection Challenge is presented. Good agreement is found with previous kinetic simulations using 
particle in cell (PIC) codes, confirming both the PIC and the Vlasov code. In the latter the complete 
distribution functions ft (k = i, e) are discretised on a numerical grid in phase space. In contrast to 
PIC simulations, the Vlasov code does not suffer from numerical noise and allows a more detailed 
investigation of the distribution functions. The role of the different contributions of Ohm's law are 
compared by calculating each of the terms from the moments of the ft- The important role of the 
off-diagonal elements of the electron pressure tensor could be confirmed. The inductive electric 
field at the X-Line is found to be dominated by the non-gyrotropic electron pressure, while the 
bulk electron inertia is of minor importance. Detailed analysis of the electron distribution function 
within the diffusion region reveals the kinetic origin of the non-gyrotropic terms. 

PACS numbers: 02.70.-c 52.25.Dg 52.65.Ff 52.25.Xz 



I. INTRODUCTION 

Magnetic reconnection is the fundamental process 
which allows magnetized plasmas to convert the energy 
stored in the field lines into kinetic energy of the plasma. 
It plays an important role in the dynamics of space and 
laboratory plasmas. In the magnetopause it allows par- 
ticles from the solar wind to enter the magnetosphere. 
Also it is believed to be the main source of energy for 
solar flares and coronal mass ejections. 

In ideal magnetohydrodynamics (MHD) the frozen-in 
flux condition prohibits the magnetic field topology to 
change. Thus reconnection depends on a mechanism that 
breaks the frozen-in condition. This non-ideal mech- 
anism is responsible for the dynamics of the diffusion 
region, where the topology change takes place. In re- 
sistive MHD with low values of resistivity an elongated 
Sweet-Parker current sheet develops, which limits the re- 
connection ratei2i£ Higher values for the resistivity, on 
the other hand, are extremely unrealistic for the astro- 
physical reconnection processes. 

In the last years it has become apparent that the 
minimal model to adequately describe collisionless re- 
connection is Hall-MHDi 4 i 5 i 6 i 7 In Hall-MHD the ion in- 
ertial length Ai is introduced as a characteristic scale 
length such that for distances smaller than Aj the ion 
dynamics can decouple from the magnetic field. Within 
the framework of the Geospace Environmental Modeling 
(GEM) Magnetic Reconnection Challenge a two dimen- 
sional magnetic reconnection setup based on the Harris 
equilibrium was studied. Simulations were done using 
MHD models with constant or localized resistivity 8-9 , 
Hall-MHD models2i2iiS, hybrid models where electrons 
are treated as fluid and ions are treated kineticallj*iiii 
and with fully kinetic particle codes ^l^H The models 
which included the Hall effect, either directly or implic- 
itly in the kinetic treatment, showed very similar results 
in terms of the global reconnection rate. The small scale 
structures of the dissipation region, on the other hand, 
varied substantially between the models. In Ref. it 



was concluded that the large scale ion dynamics control 
the rate of reconnection and that the dispersive charac- 
ter of Whistler modes is essential in the understanding 
of the results. 

These results emphasize the importance of the Hall 
term for the reconnection. At the neutral X-line, the 
magnetic field strength is exactly zero and the Hall 
term vanishes. Only the electron bulk inertia and the 
non-gyrotropic part of the electron pressure tensor can 
yield a contribution at the X-line. Kuznetsova et a/Ji 
pointed out, that the bulk inertia of the electrons can 
only contribute if the current sheet develops on elec- 
tron timescales. Vasyliunasi& was the first to emphasize 
the important role of the off-diagonal components of the 
pressure tensor in the diffusion region. Later their role 
was inv estig ated by a number of authors (see e.g Refs. 
Illlll5lll7lllal . The non-gyrotropic pressure can be un- 
derstood as originating from the meandering motion of 
the electrons, which bounce in the magnetic field reversal 
region, also known as Speiser motion^ 

In this study, we use a Vlasov code to investigate the 
GEM reconnection. In contrast to Particle-In-Cell (PIC) 
codes, Vlasov codes discretise the complete distribution 
functions of electrons and ions on a grid in phase space. 
This has the advantage of allowing more detailed inves- 
tigations of the distribution functions, since numerical 
noise is completely absent from the Vlasov-approach. 
The main focus of our investigation is the electron dy- 
namics. These can be understood in two ways. On one 
hand, one can take a microscopic point of view by looking 
at the detailed structure of the distribution function. On 
the other hand one can look at the moments of the dis- 
tribution function and their importance for Ohm's law. 
The Vlasov-approach allows us to do both with very high 
accuracy. 

In the next section we briefly present the methods used 
in our investigation. The GEM reconnection setup in- 
cluding the initial conditions and the boundary condi- 
tions is described in section ITTT1 In section HVI we discuss 
the results. Separate subsections are dedicated to the dis- 
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cussion of the contributions of the terms in Ohm's law, 
especially the off-diagonal components of the pressure 
tensor, and to the detailed discussion of the electron dis- 
tribution function. Section [V] will give a summary and 
present conclusions. 



II. METHODS 

The basis of the kinetic description of magnetic recon- 
nection are the distribution functions /fe(x, v, t) of species 
k, where k = i, e denotes ions or electrons. The tempo- 
ral behavior of the distribution function is given by the 
Vlasov equation 



^+v.V/ fe + ^(E- 
at ink 



v x B) • V v / fc = 0. 



Here qk and nik are the charge and the mass of the par- 
ticles of species k. The Vlasov equation describes the 
incompressible flow of the species phase space densities 
under the influence of the electromagnetic fields. 

The fields are solved using the Darwin approximation 
(see, for example Refs. I20l2l|) 
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where B is the magnetic and El and Et are the lon- 
gitudinal and the transverse part of the electric field, 
V x E L = 0, V • E T = 0, E = E L + E T . The Dar- 
win approximation eliminates the fast electromagnetic 
vacuum modes while all other plasma modes are still de- 
scribed. Instead of neglecting the displacement current 
completely, only its transverse part is dropped while its 
longitudinal part is kept. The elimination of the vac- 
uum modes allows larger time-steps in the simulation 
since only the slower non-relativistic waves have to be 
resolved. 

The charge density p and the current density j are 
given by the moments of the distribution function, 



P = ^2<lk J / fc (x,v)d 3 i; , 

k 

J = 5Z 9fe / vfk{*,v)d 3 v 



(5) 
(6) 



These moments couple the electromagnetic fields back to 
the distribution function. In this way the Vlasov-Darwin 
system constitutes a non-linear system of partial integro- 
differential equations. 

For the simulations we use a 2^ -dimensional Vlasov- 



we restrict the simulations to 2 dimensions in space but 
include all three velocity dimensions. In contrast to PIC 
codes, the distribution function is integrated in time di- 
rectly on a numerical grid in phase space. The integra- 
tion scheme is based on a flux conservative and positive 
scheme 2 ^, which obeys the maximum principle and suffers 
from relatively little numerical diffusion. The main idea 
of the scheme is to calculate the fluxes into and out of a 
grid cell by interpolating the primitive of the distribution 
function. The one dimensional scheme has been general- 
ized to two spatial and three velocity dimensions using 
the backsubstitution method described in Ref. |23.. The 
backsubstitution method has been shown to be slightly 
superior and, more importantly, much faster than the 
straightforward time-splitting scheme. 

The Maxwell equations in the Darwin approximation 
can be recast into a form which does not include any 
time derivatives (see, for example, Refs. I20ll2lfl . This 
makes it possible to express the electromagnetic fields 
E(to), B(to) at any time t by the density p(to) and the 
current density j (to) at the time to. No time integration 
of the fields is necessary and the conditions V • B = 
and V • E = p are always met by construction. 



III. GEM RECONNECTION SETUP 

The reconnection setup is identical to the param- 
eters of the GEM magnetic reconnection challenge. 14 
The initial conditions are based on the Harris sheet 
equilibrium 2 * in the x,y-plane 

B(y) = B tanh(|)x. (7) 

The particles have a shifted Maxwellian distribution 

(vl + (V y - V 0i , e ) 2 + vf) 



/oi, e (y,v) = n (y)exp 



2T, 



( 8 ) 

Here T i e are the constant electron and ion temperatures 
and Voi, e are the constant electron and ion drift velocities. 
The density distribution is then given by 



n (y) = n sech 2 



The Harris equilibrium demands that 
B 2 

-2- =n (T 4 + T e ) 

Zp 

, 2 Ti + T P 
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and 



Voe 

V 0i 
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(10) 
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code described in Ref. tLLl 



Here l\ -dimensional means, 



In addition a uniform background density — 0.2no 
with the same temperature T^e but without a directed 
velocity component is included. 
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The total GEM system size is L x = 25. 6A; in x- 
direction and L y = 12.8Aj in y-direction, where A^ is 
the ion inertial length A.; = c/cu p i with the ion plasma fre- 
quency is defined using the given by u> P i = {ne 2 / eQirrii) 1 / 2 . 
Because of the symmetry constraints we simulate only 
one quarter of the total system size: < x < L x /2 and 
< y < Ly/2. The sheet half thickness is chosen to be 
A = 0.5Ai. The temperature ratio is T e /Ti — 0.2 and 
a reduced mass ratio of TO,/m e = 25 is used. This re- 
duced mass ratio is consistent with the GEM setup and 
was chosen as a compromise between a good separation 
of scales (both in space and time) and the computational 
effort. The numerical expense was especially high for 
the kinetic codesAi2ii£ The reduced mass ratio was also 
used in the hybrid code a 4 ^ 1 for the sake of comparison. 
A higher mass ratio (i.e. a smaller electron mass) re- 
sults in a smaller electron skin depth which has to be 
resolved by the grid. At the same time the numerical 
time step has to be reduced to resolve the electron Lar- 
mor frequency. Both effects together imply that, in a 
spatially two-dimensional simulation, the computational 
cost increases with the square of the mass ratio for all ki- 
netic simulations independent whether thay use PIC or 
Vlasov-methods. Simulations with substantially higher 
mass ratios using the generally more demanding Vlasov 
approach, therefore, currently seem out of reach. The 
simulation is performed on 256 x 128 grid points in space 
for the quarter simulation box. This corresponds to a 
resolution of 512 x 256. This implies a grid spacing of 
Ax = Ay = 0.05Ai. The resolution in the velocity space 
was chosen to be 20 x 20 x 40 grid points. In v z direc- 
tion the grid was extended to account for the electron 
acceleration in the diffusion region. The simulation was 
performed on a 32 processor Opteron cluster and took 
approximately 150 hours to complete. 

An initial perturbation 

ij>(x, y) = ipo cos(2ttx/L x ) cos(ny j r L y ) (13) 

is added to the magnetic vector potential component 
A z . The magnitude of the perturbation is chosen to be 
tpo = 0.1£>o/Aj. The rather high values of the initial per- 
turbation generates a single large magnetic island. The 
initial linear growth of the tearing mode is bypassed and 
the system is placed directly in the nonlinear regime. The 
reason for this relatively strong perturbation is that the 
initial growth of the instabilities depends strongly on the 
electron model. In contrast to this, the GEM challenge 
demonstrated convincingly, that the later nonlinear stage 
is not sensitive to the details of the underlying model. 



IV. SIMULATION RESULTS 

As a measure of the reconnected magnetic field we use 
the difference A "J? of the magnetic vector potential com- 
ponent A z between the X-point and the O-point. Figure 
n shows the reconnected magnetic flux as a function of 
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FIG. 1: Time evolution of the reconnected magnetic flux 
throughout the simulation run. 



time. The evolution of AW agrees very well with the re- 
sults of other simulations of the GEM challenged After 
an initial small increase, the flux starts to rapidly in- 
crease at tVli rs 15. A value of AW = Bo/Xi is reached 
at t£li w 17.7. This is slightly later to what is ob- 
served in the PIC simulations. Pritchetlii^ reports a time 
tili w 15.7 where the same flux level is reached. While 
the level of saturation is comparable to the other GEM 
results, it is again reached slightly later than in Ref. |^ 
at time t£li 30. 

The upper panel of Figure shows the out of plane 
magnetic field B z at Q,it = 17.7. One can clearly iden- 
tify the quadrupolar structure generated by the Hall 
currents The peak value of the magnetic field |B Z | 
in the island structure is ph 0.17-Bo and is located at 
| a; | ~ \.7\i. This is a somewhat lower magnetic field 
than in Ref. and is also located slightly closer to the 
X line. In the course of the simulation both the size of 
the magnetic island and the magnitude of the peak mag- 
netic field increase. For this reason, the differences to 
the PIC simulations can be attributed to the small dif- 
ferences in the temporal behavior. This also explains the 
magnetic field towards the top and bottom boundaries 
of the simulation box. These fields, which can be seen 
in Figure [21 have essentially disappeared at f2,i = 20, 
approximately the same time that the peak magnetic re- 
connection rate is seen. 

The middle and lower panel of Figure [21 show the elec- 
tron and ion out of plane current densities j e%z and ji tZ at 
Slit = 17.7. While the ion current density almost exactly 
follows the ion number density rii, the electron current 
density is strongly enhanced near the X-line where the 
number densities are depleted. The thickness of the elec- 
tron current layer is smaller than the ion skin depth but 
larger than the electron skin depth. The size is deter- 
mined by the meandering motion of the electrons around 
the neutral line. In addition to this current sheet, we can 
observe thin current layers emanating from the X-line 
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FIG. 2: The out of plane magnetic field B z (upper panel), the electron out of plane current j e , z (middle panel) and the ion out 
of plane current ji tZ (lower panel) at time Qit = 17.7 



which run along the separatrix. These structures have 
not been reported in previous PIC simulations but can 
be seen in Hall-MHD simulations. 4 We point out again 
that the mass ratio in the GEM simulations was fixed 
at rrii/m e = 25 for both PIC 4 * 1 ^* 1 ^ and hybrid 4 * 1 ^ sim- 
ulations. The difference between the Vlasov code and 
the PIC code, therefore, can only be explained by the 
differences in the numerical approach. 

The electron and ion bulk velocity profiles v x (x) along 
the direction of the current sheet are shown in Figure 
|3J One can see that the electrons are ejected away from 



the X-line at super-Alfvenic speeds. These velocity pro- 
files are almost identical to those reported in Ref. [T3- In 
contrast to those results, no averaging over a finite time 
period had to be carried out here because the Vlasov 
simulations do not suffer from artificial numerical noise. 



A. Ohm's Law 

Within the GEM reconnection challenge it has be- 
come clear that the Hall-MHD model is a minimal model 
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FIG. 3: Velocity profile at fi^t = 18.1 as a function of a; at 
the location of the current sheet z — for electrons and ions. 



to understand collisionless reconnection? 4 In Hall-MHD 
Ohm's law has the form 

m dj „ „ 1 . „ 1 „ 

— -£ = E + v ; x B j x B H V • P 

ne z dt ne ne 

where the resistivity has been neglected. This is the ex- 
act electron momentum equation which can be derived 
from kinetic theory of a collisionless plasma without any 
approximations. At large scale lengths only the MHD 
terms play a role, while the Hall term and the electron 
pressure gradient can be neglected. To investigate the 
regions in which the terms of the generalized Ohm's law 
become important, we calculated the different contribu- 
tions in the whole reconnection region. 

The top panel of Figure 0] shows the inductive electric 
field E z . This inductive field is necessary for reconnection 
to take place. The region of enhanced E z is situated in a 
relatively large region around the X-line. The peak elec- 
tric field is located in an elongated area extending about 
two ion inertial lengths in the y-direction and about 4 
ion inertial lengths in the x-direction. In contrast, the 
region where the field exceeds half its peak value is almost 
circular with a diameter of about 5 ion inertial lengths. 
The middle panel of Figure 0] shows the ^-component of 
the Hall term (j x B) z . In contrast to the inductive elec- 
tric field, this quantity shows a more detailed structure. 
Two strong peaks are found left and right of the X-line. 
The peak values slightly exceed the maximum value of 
E z . These peaks coincide with the maxima of the elec- 
tron outflow velocity (see Figure |3J) . This shows that 
the Hall term is most important in the outflow regions 
where the electrons are accelerated to super Alfvenic ve- 
locities. In addition, two weaker peaks are located above 
and below the X-line, where the electrons are acceler- 
ated towards the X-line and the electron velocity starts 
to diverge from the E x B-velocity. Due to the sym- 
metry conditions, the Hall-term is exactly zero at the 
X line itself. In addition to the structure around the X- 



line, we also observe sheets of negative valued Hall-term 
along the separatrix. We attribute this to the current 
loop that generates the quadrupolar magnetic field B z . 
Away from the X-line, the electrons responsible for the 
current have to cross the separatrix back into the up- 
stream region in order to close the loop. Therefore, the 
Hall-term will have negative values along the separatrix. 
The magnitude of the Hall-term here is almost half the 
peak magnitude in the X-line region. 

The bottom panel of Figure0]shows the distribution of 
— (vi x B) . This term becomes non-zero when ions can 
move across the magnetic field lines in a region of a few 
ion inertial lengths around the X-line. Again two peaks 
can be observed in the outflow region. The peak values 
are, however, less than half of the inductive electric field. 
A striking feature in this picture is the almost circular 
ring around the X-line, where the ions become demagne- 
tized. The sheets of enhanced value along the separatrix 
are narrower than those observed from the Hall-term. 
They have the same sign as the peaks near the X-line 
and therefore partially cancel the Hall-term. 

FiguresEland|n|display the components of the electron 
pressure tensor. Although only the two mixed elements 
P xz and P yz play a role in the z-component of Ohm's 
law, the other elements are shown for completeness. The 
upper panel of Figure [S] shows the diagonal terms of the 
pressure tensor. Here we see the differences in the heating 
of the electrons in the three directions. In all three cases 
the maximum is reached in the outflow region, where the 
electron velocities are super Alfvenic. The heating in 
the x-direction, which is mainly parallel to the magnetic 
field lines, is strongest. P xx is increased mostly in the 
outflow region, with only a slight increase in the diffusion 
region. In the outflow region P yy and P zz are comparable, 
because these two directions are roughly perpendicular to 
the magnetic field lines. Within the diffusion region P yy 
and P zz are, however, different. While P yy is enhanced in 
a narrowing X type region, P zz shows a bar like structure. 
The reason for this structure of P zz may be seen in the 
acceleration of a part of the electron population in the 
^-direction. This will become more apparent in the next 
section, where the electron distributions are investigated 
in detail. 

The off-diagonal elements of the pressure tensor are 
shown in Figure [S] The magnitude of these is roughly 
one order of magnitude smaller than the diagonal el- 
ements which agrees remarkably well with the results 
of Kuznetsova et al. 11 The P xy component (top panel) 
shows a quadrupolar structure similar to the out out 
plane magnetic field B z . The P xz component has two 
extrema left and right of the X-line at the edges of the 
diffusion region in agreement with Ref. IllL In addition, 
we find two more extrema along the y = line in the elec- 
tron acceleration region and also enhanced values around 
the separatrix. Finally, P yz shows a double bar structure 
in the diffusion region and also extrema are found in the 
electron acceleration region. 

The terms of the pressure tensor that contribute at the 
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FIG. 4: 2-component of the inductive electric field E z 
2-component of v; x B (lower panel) at time Qit = 17.7 



(upper panel), Hall Term j x B (middle panel) and the negative 



X-line x = y = are given by dP xz /dx and dP yz /dy. In 
Figure[7|we have plotted the inductive electric field E z at 
the X-line over time together with the two gradients of 
the off-diagonal elements of the pressure tensor. We can 
clearly observe that the two contributions dP xz /dx and 
dPy Z /dy are roughly equal. The sum of the two shows 
are remarkable agreement with the electric field over the 
whole time of the simulation. This indicates that, at the 
X line, the bulk inertia plays only a minor role. The 
bulk inertia scales like X e /L (see Ref. ITt|) . where A e is 
the electron inertial length and L is a typical gradient 



scale length. Around the zero line of the magnetic field, 
the scale lengths of the electron dynamics are given not 
by the electron inertial lengths, but by the larger scale of 
the meandering electron motionii^iiSiSi For this reason, 
the contribution of the non-gyrotropic pressure exceeds 
the bulk electron inertia. For more realistic mass ratios 
we expect the electron bulk inertia to be completely neg- 
ligible. The pressure terms should, on the other hand, 
remain important also for higher mass ratios. Note that 
the importance of the non-gyrotropic pressure has been 
shown only close to the X-line. Away from the X-line, 




FIG. 5: The diagonal components of the pressure tensor at 
time Qit = 17.7 
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FIG. 6: The off-diagonal components of the pressure tensor 
at time Qit = 17.7 



but still inside the diffusion region, we find the Hall-term 
j x B to play a dominant role. 



To analyze the kinetic mechanism that is responsible 
for the generation of the non-gyrotropic electron pres- 
sure, we have to look at the structure of the electron 
distribution functions in the vicinity of the X-line. Fig- 
ure OH shows isosurface plots of / e (x, v) in velocity space 
for various fixed positions (x, y) in the simulation box. 
The sequence of sample points are chosen to follow a 
path from inflow region, just outside the diffusion region, 
through the diffusion region and to a point in the outflow 
region. The first panel shows the distribution function 
at (x,y) = (0.025,0.675). This is close to the symme- 
try axis x — and in the inflow region, just outside the 
diffusion region. Note that the points could not be cho- 
sen to lie on the lines of symmetry (x = or y = 0) 
since these were between the grid points in configuration 
space. Here the distribution function is slightly elongated 
in the ^-direction indicating a slightly larger pressure 
Pxx > Pyy,Pzz at this point. This moderate increase in 
P xx can also be observed in Figure |SJ While P yy and 
P zz exhibit a sharp increase at the inflow edge of the 
diffusion region, P xx rises more gradually. The distribu- 
tion at this point is totally made up of the background 
electron population. The point (x,y) — (0.075,0.275) 
lies just inside the diffusion region. In addition to the 
original background population, here two new electron 
populations suddenly appear. Both of these have a neg- 
ative w 2 -velocity. The first population is elongated in 
v x and v y direction. It is slightly bent which indicates 
a gyro motion. We interpret this population as being 
the background electrons originating from the other side 
of the current sheet. These electrons move through the 
magnetic zero line and show a bunched gyro motion when 
they enter the region on the other side with the oppositely 
directed magnetic field. Finally the second new electron 
population shows a sharp distribution with a strong nega- 
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(x = 0.025, y = 0.675) 



(x = 0.075, y = 0.275) 
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FIG. 8: Isosurface of the electron distribution function in velocity space for different positions in the simulation box. The red 
plane is v z =const, the blue plane is v x =const and the green plane is v y =const. The isosurface is drawn at f e (v x ,v y ,v z ) = 



fmax/2 where f ma x = fmax(x,y) if the maximum value of the distribution function at the position (x,y). 
ranges from Vk =-3.67 to 3.67 for all three velocity components. 



The velocity box 



tive Wz-velocity. These are the electrons which have been 
accelerated in the diffusion region. 

In the left panel on the middle row of Figure [S] the 
sample point (x,y) — (0.125,0.175) has moved closer to 
the y = line of symmetry but away from the x = 
line of symmetry. Here the populations already seen in 
the last panel become more pronounced. The gyration of 



the electrons from the opposite side of the current is now 
more apparent. This population has a distribution which 
exhibits a banana like shape, bent around the i; x -axis. 
This indicates a bunched gyro-motion. In addition, the 
distribution has been stretched to even higher negative 
v z — velocities. An slight asymmetry can be seen in the 
v^-direction due to the onsctting acceleration towards 
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the outflow region. As one moves away from the x = 
line of symmetry towards larger values of x, one must 
distinguish between those electrons flowing into the dif- 
fusion region from the ±y-directions and those electrons 
that have entered the diffusion region closer to the X- 
line. The latter population already has a strong directed 
velocity in the outflow direction. This can be seen in the 
right panel on the middle row which displays the distri- 
bution at (x,y) — (0.475,0.025). This position is almost 
on the y = line of symmetry. Again the points could 
not be chosen to exactly lie on the line of symmetry be- 
cause of the choice of the numerical grid. In this panel, 
the two populations from the ±y-directions can be iden- 
tified as two elongated blobs lying parallel to each other. 
The fact that these two blobs are well separated indicates 
that their relative velocity is larger than the thermal ve- 
locity of the electrons. For realistic mass ratios the theral 
velocity of the electrons will increase and the structure of 
the distribution function will be smeared out. Neverthe- 
less, we expect the underlying mechanisms of meandering 
motion and acceleration in the ^-direction to remain un- 
changed. The other population, which has already spent 
more time in the diffusion region, has been accelerated 
toward the outflow region. Both, the blobs and the ac- 
celerated population appear slightly tilted around the 
v y axis. This could indicate a gyro-motion around the 
newly reconnected magnetic field in the y-direction. In 
the bottom-left panel, (x,y) = (1.125,0.025), this rota- 
tion is even more pronounced. Here the different popu- 
lations start to merge and lose their individual identity. 
Finally at (x,y) = (2.425,0.025) (bottom-right panel) 
most of the structure has been lost. The temperature 
has risen considerably and a directed velocity in the in- 
direction is observed. 

The strongly structured electron distribution function 
is responsible for the off-diagonal terms of the electron 
pressure tensor. The structuring is due to the meander- 
ing motion of the electrons in the region where the mag- 
netic field approaches zero and changes signiSLSS The dis- 
tribution function is made up of a number of distinguish- 
able populations, which have a relative velocity, which is 
higher than the thermal velocity. 

V. SUMMARY AND CONCLUSIONS 

Two and a half dimensional Vlasov simulations were 
carried out on the GEM magnetic reconnection setup. 
Vlasov codes have the advantage, with respect to PIC 
codes, that they do not suffer from numerical noise and 
that the distribution function can be analyzed with high 
accuracy. This advantage is gained at the cost of a sub- 
stantially higher computational effort. We could repro- 
duce the results of other kinetic simulations carried out 
of the GEM setupi^ but were also able to calculate the 



terms in Ohm's law, especially the contributions from the 
electron pressure tensor. This shows that, although com- 
putationally more expensive, Vlasov-codes are a valu- 
able tool for investigating collisionless reconnection. The 
large scale structure of the magnetic field and the elec- 
tron and ion current densities agreed well with parti- 
cle simulations. 12 In addition we were able to identify 
some small scale structures of the electron current den- 
sity which showed enhanced values along the separatrix. 
While these structures were probably smeared out in PIC 
simulations due to the numerical noise, they resemble 
more the results seen in hybrid modelsAAi 

The analysis of the contributions to the inductive elec- 
tric field in Ohm's law show that, due to the evolution of 
the reconnection on ion timescales rather than electron 
timescales, the bulk inertia of the electron plays a mi- 
nor role. The effect of the bulk inertia will be decreased 
even more, if more realistic mass ratios are used. We 
could show that the Hall-term dominates at the inflow 
and the outflow edges of the diffusion region. Symmetry 
constraints, however, cause the Hall term and the v; x B- 
term to vanish at the X-line. Here we could confirm 
the importance of the non-gyrotropic pressure, which 
was previously investigated by Kuznetsova et The 
evaluation of the gradients of the off-diagonal pressure 
components throughout the whole time showed, that 
both dP xz /dx and dP yz /dy contribute roughly the same 
towards the electric field. The sum of the two contribu- 
tions could explain almost the complete inductive electric 
field at the X-line. 

The kinetic mechanisms responsible for the non- 
gyrotropic electron pressure were uncovered by investi- 
gating the electron distribution function. The mean- 
dering motion of the electrons in the region of the zero 
magnetic field is believed to be associated with the non- 
gyrotropic pressure>i2i2L2& The complex structure of the 
electron distribution function showed that the meander- 
ing motion is responsible. For the GEM parameters it 
is, however, not caused by the thermal electron motion 
but by the directed velocity gained in the inflow. This 
may, of course, be the result of the relatively high elec- 
tron mass. In simulations with more realistic mass ratios 
we expect the inflow velocity to be considerably smaller 
than the thermal velocity of the electrons. 
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